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CN ' An exact closed form solution for the return probability of a random walk 



on the Bethe lattice is given. The long-time asymptotic form confirms a pre- 
viously known expression. It is however shown that this exact result reduces 
to the proper expression when the Bethe lattice degenerates on a line, unlike 
the asymptotic result which is singular. This is shown to be an artefact of 

the asymptotic expansion. The density of states is also calculated. 
05.40+j 



> 






T3 

C 
O 

o 



3 



Typeset using REVTgX 



Beside being an interesting type of graph per se, the Bethe lattice (BL) is also reckoned 
as paradigm of a lattice in the limit of high dimensionality. A BL (see Fig.|l]) is usually 
defined as a set of sites connected by bonds, such that each site has the same coordination 
number and there are no closed loops. It differs from the so called Cayley tree on the fact 
that the complication arising by the boundary conditions is neglected [|TJ. 

The problem of the random walk on a BL it is not new [fjf|]. However up to now 
only asymptotic expressions were given. One of the surprising features of these asymptotic 
expressions was the difficulty arising in the interpretation of the result in the limit when the 
BL collapses into a line. 

Main aim of the present note is to derive an exact closed form solution which on one hand 
confirms previous asymptotic results, but on the other hand has the proper limit form when 
the BL reduces to a one- dimensional lattice, thus confirming the exactness of the asymptotic 
procedure and solving the aforementioned interpretation puzzle. Moreover we will provide 
an alternative solution approach with respect to the previous investigations. 

Let us start from the general master equation on the lattice: 

Px,o(t + 1) = P x ,o(t) + EK/ s ,oW - w y>x P x>0 (t)] (1) 

where P x ,o(t) is the probability density of being at site x at time t having started from site 
at the initial time to = 0. The notation y(x) means that the sum is restricted to the nearest 
neighbours y of x. In the Bethe lattice case, w x>y = 1/z where z is the coordination number 
of the lattice. 

It is convenient to introduce the generating function of P x ,o{t) (Green function): 

+oo 

^,o(A) = £A^ i0 (£) (2) 

t=o 

The fact that all points belonging to the same shell are topologically equivalent, allows to 
map the solution for the Bethe lattice onto the solution of a one-dimensional lattice with a 
defect. Therefore the Green equation takes on the Bethe lattice the form: 

P ,o(A) = AA,o(A) + 1 (3a) 



P n ,o(A) = ^P„_i,o(A) + KZ - l) P n+ i,oW (3b) 

for the 0— th and n— th shell respectively. Here P n ,o{t) refers then to the probability of being 
in the n— th shell at time t having started from the seed 0. 

The solution of eq. (|3a| , |3b|) is considerably simplified by noting that the ratio 
-Pn+i,o(A)/-P n ,o(A) is independent on n due to the homogeneity of the lattice and the par- 
ticular boundary conditions [|J. It is then a simple matter to solve the quadratic equation 
coming from (J3bj) and substitute the root, which has finite value in the A — ^ limit, into 
([3a|) with the result: 

P 2(z-l)/z 

This result was previously obtained by Cassi by a different procedure. 

It is important to notice that for z = 2 this expression reduces to the well known || 
result of the generating function for the one dimensional lattice. It is also worth mentioning 
that since the critical value for the fugacity A c = z/2y/ ' z — 1 > 1 for z > 3 the generating 
function P ,o(A) is always real and finite for A < 1 and therefore a random walk on the Bethe 
lattice cannot be critical. 

Upon serie expansion of eq.(F§) and using the definition (E3), one gets after some algebra: 



P (9,1 (s-^ f V^V y (2p + 2t)! z-1 

flw(a) = -—(-—) Y. g(p + t + mp + ty (—y (5) 

and Po t o(2t + 1) = for t > with -Po,o(0) = 1. After some manipulations this expression 
can be cast in the following closed form: 



W t) = ^(^ V+^rV+i) > m + 1/2 ' *• * + 2A(Z ~ 1),z2) (6) 

where T(t) is the gamma function and 2 Fi(a,P, 7, z) is the Gauss hypergeometric function 



For large t one can use the property [[?] 

2 F 1 (a,P,P,z) = (l-z)- a (7) 



valid for arbitrary (3 and the Stirling approximation for the gamma function || to find, at 
the leading order in t » 1: 

^«'~ 1 3^ S/2expHln( 57^T )) (8) 

which confirms the asymptotic result derived in ref . [[| . Note that in the limit z — > 2 (when 
the Bethe lattice degenerates on a line) the expected asymptotic behaviour Po,o(t) ~ t^ 1 ^ 2 
is not recovered both because the prefactor becomes singular and because the (universal) 
power law is not the correct one. 

We are now in the position to show that this is an artefact of the asymptotic expansion 
stemming from the fact that the limit t — > oo and z — * 2 do not commute. Indeed if we set 
z = 2 in the exact result (|6|) and use the fact that |J: 

. 2 F l(a ,ft 7 ,l) = 'S 7)r(7 - a ~«! (9) 

r(7 - a)r(7 - p) 

we get: 

p °»v = "tf ^+TAti) (io) 

which is the well-known result of the one-dimensional case ||. This latter result could in 
fact be derived by starting from eq.@ for z = 2 and using a procedure similar to the one 
employed here. 

It is also possible to compute the density of states, associated to the master equation ([[]), 
on the BL. This was not previously calculated. Indeed using the transformation procedure 
between discrete and continuum times described in ref [0], it is not hard to see that the 
the Laplace transform of the return probability is related to the generating function P 00 (A) 
through: 

Po,oH=APo,o(A)U=i/(i+u,) (11) 

where ui is the Laplace variable conjugate of time. Then, using eq.(f|) and ([LTD, one obtains: 

A,M = r^ 2( r 1)A " < 12 > 



The density of state p(e) is then well known || to be given by the analytical continuation: 

p (e) = --lmP ofi (-e + iO + ) (13) 

TV 



Then, in the present case, the result of the analytical continuation is: 

if 2e - e 2 - (z - 2) 2 /z 2 > 



P(e) 



f ^ v^MpVZ if 2e _ e 2 _ {z _ 2)2/^2 

2tt 2e- e 2 V ) I ^ 



otherwise 



which for z = 2 reduces to the well known formula of the density of states of a one- 
dimensional lattice |1(J . FigfJ shows the comparison between the case of the one-dimensional 



lattice (z = 2) and the BL (z = 3,4). 

Summarizing we presented an exact closed form solution for the return probability and 
the density of states on the BL. Previous asymptotic results were confirmed and explained 
in terms of this solution. Although other relevant quantities, beside the ones presented here, 
could, in principle, be obtained, the required algebra becomes rapidly very involved. This is 
beyond the purposes of the present work whose main objective was the confirmation of the 
asymptotic results, along with the removal of the inconsistency contained in them. As a by- 
product of our investigation we also presented an alternative simplified solution procedure 
with respect to the previous two approaches. 
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FIGURES 
FIG. 1. Example of a Bethe lattice with coordination number z = 3. Sites indicated with the 

same numbers n = 0, 1, 2, ... belong to the same shell. 

FIG. 2. Density of states associated to the master equation in the case of a one-dimensional 
lattice (z = 2) and of a Bethe lattice (z = 3,4). 



